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Abstract. We use a master equation approach based on the Poland- Scheraga free 
energy for DNA denaturation to investigate the (un)zipping dynamics of a denaturation 
wedge in a stretch of DNA, that is clamped at one end. In particular, we quantify the 
blinking dynamics of a fluorophore-quencher pair mounted within the denaturation 
wedge. We also study the behavioural changes in the presence of proteins, that 
selectively bind to single-stranded DNA. We show that such a setup could be well- 
suited as an easy-to-implement nanodevice for sensing environmental conditions in 
small volumes. 

PACS numbers: 87.15.-v, 05.40.-a, 82.37.-j, 87.14.Gg 
1. Introduction 

During the last decade or so, technical progress in detection and manipulation of 
single molecules and their dynamics has snowballed. By fluorescence spectroscopy, 
optical tweezers, atomic force microscopy, or patch clamp techniques, for instance, it is 
possible on the single molecule level to probe the opening and closing dynamics of local 
denaturation zones in a DNA molecule 1], to study the binding of single-stranded DNA 
binding proteins to overstretched DNA 00], to stepwise disrupt domains in a protein [I], 
or to monitor the passage of a biopolymer through a nanopore in a membrane [3 El, just 
to name a few. This experimental progress is accompanied by advances in the theoretical 
understanding of fundamental processes relevant on small scales such as the Jarzinsky 
relation connecting measurements of the nonequilibrium work needed, e.g., to stretch 
an RNA segment [7j, to the difference in the corresponding thermodynamic potential 
[H]; or the entropy production along single trajectories exposed to stochastic forces 0. 
Given the novel possibilities to synthesize supramolecules with topologically confined 
mechanical units JU], entropy-based designer molecules such as molecular muscles were 
proposed [HUE], and new possibilities to produce dynamic nanosensors discussed [T3] . 

In what follows, we explore the dynamics of a stretch of double-stranded DNA 
(dsDNA), that is clamped at one end but allowed to denature from the other, as 
sketched in Figure [U Internal bubble formation, in comparison, is suppressed by a 




Figure 1. Schematic of the end-denaturing of a double-stranded DNA molecule, that 
is clamped at one end; here, by attachment to a wall. The number of denatured base- 
pairs is m, and the overall length of the DNA segment is M . Selectively single-stranded 
DNA binding proteins from the surrounding solution bind to the denatured portion 
of the DNA. Once bound, the SSBs prevent closing of the denaturation wedge. In 
reality, the clamping can be achieved by sealing the denaturation zone created by AT 
base-pairs with a stretch of more stable GC base-pairs, compare reference pp. 



Boltzmann factor cr ~ 10~ 5 . . . 10~ 3 |14j . and this effect can therefore be neglected, 
see below. In the setup we have in mind an individual base-pair is tagged with a 
fluorophore-quencher pair, see Figure El Once separated, the dye starts to fluoresce, 
and the resulting blinking can be detected [T5j . Our aim is to establish a quantitative 
description of such a molecular beacon. We will use a master equation approach to 
describe the sequential unzipping and zipping dynamics of the denaturation wedge, also 
taking into account the possible presence of proteins, that selectively bind to single- 
stranded DNA (ssDNA). These single-stranded DNA binding proteins (SSBs) occur in 
the cells of most organisms. Here, we show that the autocorrelation function for the 
blinking depends on the concentration of SSBs and the binding strength between the 
single strand and the SSBs (the latter varies, for instance, with salt concentration) as 
well as temperature. The denaturation wedge depicted in figure El therefore acts as a 
nanosensor that can be probed on the single molecule level using fluorescence techniques. 
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Figure 2. Schematic of the molecular beacon setup. A fluorophore starts to fluoresce 
in the incident laser light, once the denaturation wedge moves the fluorophore apart 
from the quencher. 

2. Experimental setup: Quantifying the blinking dynamics of the molecular 
beacon 

Before setting out to describe our general theoretical scheme, we first describe the 
experimental setup we have in mind in some detail. We consider modelling the blinking 
behaviour of a fluorophore-quencher pair mounted on the denaturation wedge as shown 
in Figure El This setup, similar to the ones described in References [TJ [T3] works as 
follows. As long as the dsDNA is intact, fluorophore and quencher are in close proximity. 
Once they come apart from one another when the denaturation wedge opens up, the 
incident laser light causes fluorescence of the dye. The on/off blinking of this "molecular 
beacon" can be monitored in the focus of a confocal microscope. 

The blinking renders immediate information about the state of the base-pair, that 
is tagged by the dye-quencher pair. Blinking, that is, indicates that the base-pair is 
currently broken. It is therefore advantageous to define the random variable I(t) with 
the property 

^, , I if base-pair at x = xt is closed ,^ 
| 1 if base-pair at x = xt is open ' 

and we are interested in the behaviour of the autocorrelation function 



A(t) = (J(t)/(0))-(J) eq2 , 



(2) 
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where (I)eq is the (ensemble) equilibrium value. Given the fact that the formation 
of an internal denaturation bubble is connected with a rather high initiation barrier 
o"o — 1CT 5 . . . 1CT 3 (corresponding to 7. . . 12 ksT at room temperature), such bubbles 
are much less probable than denaturation from the undamped end, we focus only on the 
end-denaturation. Therefore, a base-pair at x = xt is open if m > xt, see figure El A 
word on the interpretation of the average (I(t)I(0)) is in order. Denoting by p(I, t; I', 0) 
the associated joint probability density that the tagged base-pair is in state / at time 
t given that it was in state I 1 at initial time t = 0, we can rewrite the autocorrelation 
function as 

(i(t)/(o)> = j p( j > *; J, > °) J/ = *; !> °)- ( 3 ) 

This is but the survival probability density for a denaturation wedge, i.e., the probability 
density that the base-pair is open at time t, given that is was open at system preparation 
at time t = 0. 

In the remaining part of this study we present a general scheme for calculating 
the opening-and closing dynamics for the physical system presented in figures Q and El 
The focus is on presenting a scheme which allows for the calculation of the measurable 
quantity A(t) defined above. We note, however, that our scheme is sufficiently general 
to allow for straightforward derivation of other experimental quantities, such as the 
binding dynamics of SSBs in optical tweezers overstretching setups [21 El QUI- m that 
case, the Boltzmann factor for opening a base-pair, u = exp(/3e) (e being the binding 
energy of a base-pair, see below) becomes modified to u = exp (j3e — T9 ), where T 
is the external torque exerted by the optical tweezers setup, and 9 = 27r/10.35 is the 
relaxed dsDNA-twist per base-pair. 



3. Master equation for end-unzipping 

Denote by P(m,n,t) the probability distribution that there are m broken base-pairs 
and n bound SSBs at time t. As m and n are the slow variables of the system, their 
dynamics can be described in terms of the master equation 
oPim Ti t] 

v ' l > = t+(m - 1, n)P(m - 1, n, t) + t"(m + 1, n, t)P(m + 1, n, t) 

dt 

— [t + (m, n) + t~(m, n)\ P(m, n, t) 

+ r + (m, n — l)P(m, n — 1) + r~(m, n + l)P(m, n + l,t) 

— [r + (m, n) + r~(m, n)] P(m, n, t), (4) 

compare the discussion in References [TTJ HHJ HH|. Mediated by the transfer rates t 1 * 1 
and r ± , the population of a given state (m,n) is continuously changed by (un)zipping 
a further base-pair and (un)binding of an SSB. To complete the master equation (JU), 
we need to specify the boundary conditions: The clamping to the right at base-pair 
x = M + 1 ensures that no further unzipping beyond base-pair M occurs, namely, 

t+(M,n) = 0. (5) 



m = ( n + 1 ^ X ,n = n max (m) - 1 ) = 0. (8) 
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Moreover, if both branches of the denaturation wedge are fully occupied by SSBs, i.e., 
the maximum number of SSBs, 

Tfl ' 

n mH (m) = 2 — , (6) 
. A J 

is bound, then the base-pair at the zipping fork is not allowed to close: 

m = ^,n = n max (m)^) = 0. (7) 

Also, if only one of the two branches of the wedge is fully occupied the zipper cannot 
close: 

t I m, = - — 

2 

Similar to the (un)zipping rates t 1 * 1 , we impose the boundary condition 

r+ (m,n max (m)) = 0, (9) 

i.e., once the denaturation wedge is completely occupied, no additional SSB is allowed 
to bind; and when n = SSBs are bound, no further SSB can detach: 

r -(m,0) = 0, (10) 

The configuration lattice showing the allowed moves is illustrated in Figure H3 Empty 
arrow heads indicate forbidden moves. 

The general solution of the master equation (pjj can be obtained through the ansatz 

P(m, n,t) = Y^ c p Qp(m, n)e~** , (11) 
p 

corresponding to an expansion in eigenmodes. Here, the expansion coefficients c p 
are determined by the initial condition. Inserting above expansion into equation (jjj) 
produces the eigenvalue equation 

— r]pQ p (m, n) = t + (m — 1, n)Q p (m — 1, n) + t~(m + 1, n)Q p (m + 1, n) 

— [t + (m, n) + t~(m, n)l Q p (m,n) 

+ r + (m, n — l)Q p (m, n — 1) + r~(m, n + l)Q p (m, n + 1) 

— [r + (m, n) + r~(m, n)] Q p (m,n) (12) 

for the pth eigenmode. The concrete forms for the transfer rates are defined in the next 
section. 



4. Partition function and transfer rates 

In order to obtain the transfer coefficients t 1 * 1 and r ± , we first define the partition function 
3f(m, n) for the end-denaturation fork of the dsDNA and the SSBs bound to its two 
branches. To this end, we note that we can decouple the partition coefficient 

iF(m, n) = ^ DNA (m)^ SSB (m, n) (13) 
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Figure 3. Configuration lattice showing the possible configurations □ of the system 
in the (m, n)-plane. The empty arrow heads represent forbidden moves due to the 
boundary conditions. The maximum number of bound SSBs is N = 2[M/X]. 



into the contributions jg° DNA counting the degrees of freedom of the DNA molecule, and 
the contribution of the SSBs. According to the Poland-Scheraga model of DNA 

denaturation, we have 



2? 



DNA/ 



m 



u 



(14) 



Here, u = exp (j3e) is the Boltzmann factor for opening a base-pair, i.e., the activation 
needed to overcome the free energy barrier e for opening an additional base-pair, e 
combines two appreciable contributions from enthalpy cost and entropy gain on breaking 
the base-pair, that almost cancel such that e is of the order of a k^T at physiological 
temperature and salt: At 37°, u ~ 0.6 in a zone of AT base-pairs, and u ~ 0.2 in a GC- 
domain [14J. We note that the parameter u is sensitive to salt concentration. Here, we 
consider homopolymer zones of either AT or GC; for the treatment of a heteropolymer, 
see Reference [2D]- As already mentioned, we also neglect the formation of internal 
denaturation bubbles within the dsDNA stretch, that would require the crossing of the 
initiation barrier cr . This contrasts our previous studies j2B H3 CHI EOl E2j , in which we 
focused on internal bubbles, preventing end-denaturation wedges by clamping at both 
ends. Note also that in the case of the end-wedges, we do not have to take care of the 
entropy loss due to loop formation as present in internal denaturation bubbles. 
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(17) 

n /< n max/2 



The contribution from the SSBs has the form 

^ SSB (m, n) = K n Q SSB {m, n), (15) 

sharing the energetic component 

k = c K eq (16) 

per bound SSB, where Co is the SSB concentration in solution and K eq = 
voexp (\Eund\/kBT) the equilibrium binding constant of SSB-binding, with the typical 
volume v of the SSBs and the binding energy E hin d', and the number of possible ways 
of arranging n SSBs on the two arches of the denaturation fork, both being of length m 
PUDS]: 

n n-ri'<n max /2 

n SSB (m, n)=J2 ^V, "> SSB (m, n - n') 

n'=0 

This counting allows to have a different number of SSBs on the two arches. The degrees 
of freedom of distributing n SSBs each covering A bases, on an arch of size m is 

oa R . . (m — (A — l)n\ (m — [A — lln)! 

w SSB (m;n)= ' = ,/ X V • (18) 

\ n / n!(m — Anj! 
From the partition function, we can now come to specify the transfer rates. As a 
first requirement we impose that the rates obey the detailed balance condition [23J 

r + (m, n — l)3f(m, n — 1) = r~(m, n)3?(m, n) (19) 

for SSB (un)binding, and 

t + (m — 1, n)3?(m — 1, n) = t~(m, n)3f(m, n) (20) 

for base-pair (un)zipping. Detailed balance guarantees that the probability distribution 
P(m,n,t) for long times relaxes towards the Boltzmann distribution. However, the 
detailed balance condition does not uniquely specify the rates. Using the remaining 
freedom of choice [23 El, we settle for the following forms. In the case of SSB 
(un) binding, we choose 

r+ m, n) = (n + 1 )q—^f ^ = n + 1 g« V B / 21 

for the binding rate, and 

r~(m, n) = nq (22) 

for unbinding. The unbinding rate, that is, is proportional to the number of bound 
SSBs, as it should. Apart from this n-factor, we choose it to solely depend on the 
constant unbinding rate q. Conversely, the binding rate includes, apart from this rate q, 
the binding strength k, the relative number of degrees of freedom given by the ratio of 
the Q SSB factors; and the factor (n + 1). The factor (n + l)f2 SSB (m, n + l)/^ SSB (m, n) 
accounts for the combinatorial number of ways to put an additional SSB of size A 
onto either arch of the denaturation fork. In the limit A = 1 the resulting expression 
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for r + (m, n) becomes qK(m — n) and thus designates the number of free binding sites, 
compare the discussion in [T71 ■ 

Similarly, we choose a completely asymmetric form for the base-pair (un) zipping 
rates. Thus, for unzipping of an additional base-pair we use $ 

t + (m, n) = t+(m) = k = hi, (23) 

carrying the full Boltzmann factor u, apart from the rate constant k. The zipping rate 

contains all information related to the interplay with the number m of bound SSBs, 
and is proportional to the probability that the base-pair next to the denaturation 
fork is unoccupied. This choice realistically describes the fact that a region almost 
fully occupied with SSBs is less likely to decrease in size. We note that in general 
t~(m, n) < k, and that t~(m, 0) = k. Note also that we can conveniently define a 
dimensionless parameter 

7 = ! (25) 

that measures the competition between SSB (un)binding and base-pair (un)zipping. 

5. Results: Blinking autocorrelation of the molecular beacon 

Having defined the dynamics of the end-denaturation wedge in terms of the numbers 
of broken base-pairs, m, and of bound SSBs, n, we now proceed to calculate the 
autocorrelation A(t) for the blinking behaviour of a fluorophore-quencher pair mounted 
on the denaturation wedge as shown in Figure El Following the results from reference 
|2Uj . we can express A(t) according to equation (J2J) in the form 

m = J2 T r 2e ~ t/Tp i ( 26 ) 

with relaxation times r v = t]" 1 , and where 



m,n 



(27) 

It should also be noted that, in order to be detected experimentally, it might be 
necessary that the fluorophore-quencher pair should be separated further than by the 
separation provided by solely opening the base-pair x T . Thus, if A additional base- 
pairs need to be broken before fluorescence occurs, the quantity T p will be modified to 

T P = Em.n Qp( m ' U ) ^ ■ 

t We here, for simplicity, neglect any m-dependence of the rate constant k, compare the discussion 
about the 'hook'-effect in references |17l I18| . 
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Alternatively, we can rewrite the autocorrelation function A(t) according to the 
spectral decomposition 

A(t)= [°° f(r)exp(--)dr. (28) 



Jo \ T / 

This way, we obtain the weighted spectral density ("relaxation time spectrum") 

This distribution f(r) renders information about the contributions of individual 
relaxation modes to the autocorrelation function A(t), see also below. 

For the two limiting cases corresponding to absence of SSBs, and to fast (un)binding 
dynamics of the SSBs, we obtain analytical solutions for A(t) from the master equation 
(@J below. The general case is solved numerically along similar lines as developed in 
references [T7| . 

5.1. Absence of SSBs 

If no SSBs are present, the (un)zipping rates t ± assume the simpler form 

t + (m) = t+(m,0) = ku, (30) 

and 

t~(m) = t + (m, 0) = k, (31) 

corresponding the transfer rates for an asymmetric random walk. This contrasts the 
non- linear form for internal bubbles that includes the loop entropy loss ^3 The 
master equation reduces to the one-variable form 

8P{ \' Q = t + (m - l)P(m - 1, t) + t~ (m + l)P(m + 1, t) 
at 

- [t + (m) + r (m)] P(m, t) (32) 

with the eigenvalue decomposition P(m, t) = J2 P c p Q p (m)e . Thus, we obtain the 
eigenvalue equation 

- r] p Q p {m) = t + (m - l)Q p (m — 1) + t~(m + l)Q p {m + 1) 

- [t + (m)+r(m)]Q p (m), (33) 

with the obvious boundary conditions t~(0) = 0, and t + (M) = 0. The solution of the 
reduced eigenvalue equation for Q p can be obtained in terms of orthogonal polynomials 
(or Chebyshev polynomials) according to References [T7| l2^j: 

Q p {m) = {sin [(m + l)u) p ] - u~ 1/2 sin [mw p ]} , (34) 

Sill LU p 

where the eigenvalues become 

r} p = k [u + 1 — 2-u 1 ^ 2 coscjp] (35) 
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with 

u v = — — -■ (36) 
v M+ 1 v ; 

We notice that the relaxation times f p = l/f) p fulfil the inequalities 

r min (u) = k~ l (1 + u 1 ' 2 )' 2 <r p < k- 1 (1 - u x l 2 Y 2 = r max (u). (37) 

Thus, at the melting transition where u — 1, and for an infinitely long DNA segment 
M —>■ oo the longest relaxation time diverges, i.e., the denaturation wedge has a 
diverging lifetime, as one would expect. The correlation function A(t) corresponding to 
the results above are shown by the black dash-dotted curves in figure |U 

We note that conditions similar to the absence of SSBs are fulfilled if the 
concentration of SSBs is very low, k — > 0, or if — > 0, compare the existence of a 
kinetic block to SSB-binding in DNA-bubbles El EH EE! ■ 

5.2. Fast (un)binding dynamics of SSBs 

In case the (un)binding dynamics of the SSBs is much faster than the typical base- 
pair (un)zipping rates, 7 ^> 1, we can adiabatically eliminate the degrees of freedom 
corresponding to the SSB dynamics |2Sj- Ensuing are an effective free energy landscape, 
dressed by the SSB- interactions with the denaturation wedge, and a reduced one- variable 
master equation of the type (|3*2j) . but with the following (dressed) rate coefficients (T7j: 

n max (m) 

t ± M= E t + K-)^^, (38) 

where 3?{m) = 3?{m, n). When calculating the two dressed rates t , it is important 
to consider the boundary conditions ©, ©) and ©• From these we also deduce the 
modified boundary conditions t + (M) = and t~(0) = 0. The corresponding result for 
the autocorrelation function are shown by the dashed curves in figure EJ 



5.3. General case 

In the general case when the dynamics of SSB (un)binding and base-pair (un)zipping 
occur on comparable time scales, 7 ~ 1, the eigenvalue equation (JT2J) has to be solved 
numerically; for detailed elaboration on the procedure, see reference jT7]. In figure^ we 
display some typical examples for the relaxation time spectrum and the autocorrelation 
function for various values of SSB binding strength k and rate ratio 7. 



5-4- Discussion 

We first note that the autocorrelation function A(t) in the semi- log plot is non- 
exponential, corresponding to a multimodal relaxation behaviour, that was also observed 
in experiment pQ. This is reflected in the rather broad distribution of relaxation times 
shown in figure 0] In the binned data for the relaxation times, the first bin contains the 
major portion of the relaxation contributions for all cases. This describes the relaxation 
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Figure 4. Relaxation time spectrum /(r) and autocorrelation function A(t) for 
different values of the binding strength k and the ratio 7 = q/k between the unbinding 
and the unzipping rate constants, as indicated in the graphs. The dashed curves 
correspond to the approximation of adiabatical elimination (fast binding) . In the plots 
for the relaxation time spectra, /(r), the data were binned (the black dashed lines 
at the top of the graphs show the bin size). We chose u = 0.6, A = 5, Xt = 1, and 
M = 30. Notice that strongly binding SSBs increase the relaxation time by orders of 
magnitude (logarithmic abscissa in the A(t) plots). 
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Figure 5. Longest relaxation time r re i ax as a function of SSB binding strength k. 
Note the dramatic variation of r ro i ax (logarithmic ordinate). We chose u = 0.6, A = 5, 
and M — 30. A maximum occurs at around k = 1. 

of the denaturation wedge itself, i.e., the relaxation due to the base-pair (un)zipping 
process. In the plots for A(t), this corresponds to the first relaxation shoulder located 
at a few inverse zipping time units, k~ x . 

In absence of SSBs, our coarse-grained plot of the relaxation time distribution 
shows only one contribution, corresponding to only one relaxation shoulder in the graphs 
for A(t). With increasing binding strength k, a second relaxation shoulder in A(t) is 
building up. This is due to the relaxation of SSB (un)binding. Let us discuss the 
somewhat involved behaviour for the various values of k in the case A = 1, for which the 
occupation fraction of SSBs on the ssDNA branches of the denaturation wedge becomes 
/ = k/(1 + k) [1^]. We see that k — 1 leads to an occupation fraction / = 1/2. In that 
case, we would expect the relaxation time for SSB occupation to be the longest. Both 
for increasing and decreasing k, the fraction of vacancies and bound SSBs is decreasing, 
respectively, so that the exchange between both species becomes faster, and A{t) decays 
quicker than for k = 1. However, for k < 1, the number of bound SSBs is smaller 
and mostly a certain binding site is vacant, and the relaxation therefore quicker than 
for k > 1, for which a site is mostly occupied. Indeed, these trends can be observed 
in the plotted examples: The relaxation time contribution is the longest in the case 
K — 1. This observation is corroborated in figure El illustrating the behaviour of the 
longest relaxation time as function of binding strength k: the maximum close to k ~ 1 
is distinct. 

In dependence of the ratio 7 = q/k of rates q for SSB unbinding and k for base-pair 
unzipping, the increase in the relative speed, 7, of the SSB dynamics is immediately 




Blinking statistics of a molecular beacon 



13 



evident from the shift towards shorter r in the relaxation time spectrum. In the plots 
for A(t), the results in presence of SSBs approach the adiabatic approximation; in the 
case 7 = 2, virtually no difference between the full result and the fast binding limit are 
visible in A(t). 

6. Conclusions 

We investigated by means of a master equation approach in detail the dynamics 
of end-denaturation of a clamped stretch of homopolymer DNA. The transfer rates 
were obtained from the partition coefficients based on the Poland-Scheraga model of 
DNA denaturation. We chose the rates t^ 1 and ^ for base-pair (un)zipping and SSB 
(un)binding, such that detailed balance is fulfilled, i.e., thermal equilibrium reached for 
long times. We furthermore chose these rates in a fully asymmetric form guaranteeing 
that base-pair unzipping is proportional to the Boltzmann factor for breaking the base- 
pair and the fundamental rate k. For base-pair zipping, the rate is given by k in absence 
of SSBs, whereas in presence of SSBs it gets dressed by the combinatorial ratio slowing 
down the zipping when SSBs are bound (essentially, by the probability that the base- 
pair next to the zipping fork is vacant), and preventing closing of the wedge once it 
is fully occupied. Similarly, we chose the (un)binding rates such that the unbinding is 
specified by the fundamental rate q times the number of bound SSBs. Binding is biased 
by the binding strength k and the probability of adding another SSB to the two arches 
of the denaturation wedge. 

Mounting a fluorophore- quencher pair at position xt on either arch of the 
associated denaturation wedge, a molecular beacon is built, whose blinking dynamics, 
corresponding to the open or closed state of the base-pair xt, is described by the random 
variable /. For this quantity, we obtained the autocorrelation function A(t) and the 
associated spectrum of relaxation times. The quantity A{t) can be measured directly 
in experiments (see, e.g., references [UHS])- The predicted behaviour of A(t) and the 
spectrum /(r) shows a multimodal behaviour with two pronounced relaxation shoulders 
in the presence of SSBs. As our model involves physical parameters that are known for 
given solvent conditions (it, k) or can be determined independently (k, q), our model 
becomes fully quantitative, and can be used to devise future experimental setups. 

The molecular beacon setup we propose represents the basis for an interesting 
nanosensor. Constructing a small stretch of AT base-pairs and clamping them with a 
few GC base-pairs at one end (compare the construction used in reference pQ), such 
a nanosensor would be of the linear size of some 20 nm. A low concentration of such 
nanosensors would therefore be sufficient to probe for presence of SSBs, salt conditions, 
or similar, in small probe volumes as, for instance, encountered in gene arrays. The 
rather large changes in the relaxation time spectrum invoked by the presence of SSBs (or, 
by varying u, corresponding to temperature or salt changes, not shown here) corroborate 
that this kind of sensor could actually be rather sensitive. We note that instead of 
the conventional fluorophores, that bleach rather quickly, longer-lived quantum dots or 
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plasmon resonant nanoparticles are now available, see, for instance, references |2b[ I27j. 

Whereas our calculations are for a homopolymer denaturation zone, the simplest 
possibility in view of designing a nanodevice, it is possible to extend the model to 
a heterogeneous sequence, see reference [20J. However, the numerical evaluation of 
the corresponding master equation may become challenging, especially, when longer 
denaturation zones and small SSBs are employed. An efficient alternative for such more 
involved cases is provided by stochastic simulations techniques such as the Gillespie 
algorithm, as studied recently in the context of denaturation fluctuations in DNA [22] • 
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